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ABSTRACT 



In a recent paper (Svanda et al., 2008, A&A 477, 285) we pointed out that, based on the tracking of Doppler features in the full-disc 
MDI Dopplergrams, the active regions display two dynamically different regimes. We speculated that this could be a manifestation 
of the sudden change in the active regions dynamics, caused by the dynamic disconnection of sunspots from their magnetic roots as 
proposed by Schiissler & Rempel (2005, A&A 441, 337). Here we investigate the dynamic behaviour of the active regions recorded 
in the high-cadence MDI data over the last solar cycle in order to confirm the predictions in the Schiissler's & Rempel's paper. We 
find that, after drastic reduction of the sample, which is done to avoid disturbing effects, a large fraction of active regions displays a 
sudden decrease in the rotation speed, which is compatible with the mechanism of the dynamic disconnection of sunspots from their 
parental magnetic structures. 
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1. Introduction 

The dynamics of active regions in the solar photosphere and in 
the close sub-photospheric layers is closely related to the so- 
lar dynamo process. The magnetic flux emerges from beneath 
the surface and forms sunspots and other active phenomena. At 
some point, the magnetic field starts to disperse, the sunspots 
disappear, and the active region decays, forming the surge-like 
structures of the trailing p olarity expanding to the solar poles 
dBumba & Howard |1965[). where they c ontribute to the solar 
field reversals (see e.g. Wang et aTlll989l) . 

The measurements of the dynamic behaviour of active re- 
gions can bring forth new insights into what is happening to 
the magnetic field under the surface. Many authors have stud- 
ied the variations in the rotation speed of sunspots. The rotation 
of the suns pots in relation to their morphological type was stud- 
ied by e.g jBalthasar et all (1 1986b . who found that more evolved 
types of sunsp ots (E, F, G, and H ty pes) rotate slower than less 
evolved types. [Ruzdiak et all d2004l) investigated the Greenwich 
Photoheliographic Results for the years 1874-1976 and found 
clear evidence for the decel eration of the sunspo ts in the photo- 
sphere with their evolution. iHerdiwijaval (|2002) found that the 
leading part of a complex sunspot group rotates about 3 % faster 
than the following part. The dependence of the rotation rate of 
sunspots on th eir size and position in th e bi-polar region was 
investigated bv uTSilva & Howard! i 19941) . Their results showed 
that smaller spots rotate faster than large ones. The authors ex- 
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plained the observed behaviour through a subtle interplay be- 
tween the forces of magnetic buoyancy and drag, coupled with 
the Coriolis force acting on rising flux tubes. 

The interpretations used in the explanation of the observed 
dynamics of active regions often involve the deep anchoring of 
the magnetic structures. From helioseismic inversions, we know 
that throughout the convective envelope, the angular velocity 
contours at mid-latitudes are nearly radial. Near the surface, at 
the top of the convection zone, there is a layer of a large radial 
shear in the angular velocity. At low and mid-latitudes, there is 
an increase in the rotation rate immediately below the photo- 
sphere that persists down to r ~ 0.95 R Q (where r is the radial 
coordinate). The angular velocity variation across this layer is 
roughly 3 % of the m ean rotation rate, and accordi ng to the he- 
lioseismic analysis of lCorbard & Thompsonl d2002l) . the angular 
velocity a> decreases within this layer approximately as r~ 1 . At 
higher latitudes, the situation is less clear. The changing depth 
of the magnetic roots, e.g., by the rising motion, could explain 
the deceleration of the sunspots during their evolution. 

One of the important properties observed in the bi-polar 
sunspot groups is their tilt with respect to the zonal direction 
(lHale et all 119191) . This tilt is believed to be generated by the 
Coriolis force, which acts to twist the flux-loop as it rises through 
the convection zone. The obvious question is, what happens to 
the tilt when the flux eruption has ceased, because this implies 
that the twisting Coriolis force no longer exists? If the magnetic 
field remains connected to the deep toroidal field at the base of 
the convection zone, then the tilt should slowly diminish through 
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the action of the magnetic tension. However, this behaviour is 
not observed, which challenges the solidity of the connection to 
the deep magnetic field. 

A number of theories h ave been suggested dFan et all 1 19941: 
ISchiissler & Rempe"ll2005l) indicating that bi-polar magnetic re- 
gions may become disconnected from their magnetic roots and 
form an isolated island-like feature. The mechanism is based 
upon the buoyant upfiow of plasma along the field lines. Such 
flows arise in the upper part of a rising flux loop during the final 
phase of its buoyant ascent towards the surface. The combina- 
tion of the pressure build-up by the upfiow and the cooling of 
the upper layers of an emerged flux tube by radiative losses at the 
surface leads to a progressive weakening of the magnetic field at 
depths of several Mm. When the field strength has become suf- 
ficiently low, convective motions ablate the flux tube into thin, 
passively advected flux fragments. The mechanism provides for 
a dynamical disconnection of the emerged part from its parental 
magnetic structure. This instant should be observed as a change 
in the dynamic regime, as the "floating island" does not reflect 
the deep dynamics any longer. We note that even in the models 
which assume that the surface magn etic activity is a consequence 
of the shallow dynamo action (e.g. ISchatterJ. l200"9h . the separa- 
tion of the surface local magnetic field from the close subsurface 
one is required, in order to allow the magnetic flux to disperse 
towards the solar poles. The model of lSchattenl d2009l) however 
does not provide the predictions for the systematic dynamics of 
active regions in the photosphere. 

In the recent paper of ISvanda et a l. (2008) we found that 
sunspots in the equatorial region seem to show two differ- 
ent dynamical regimes. In one regime, the fast regime, the 
sunspots displayed an almost constant velocity of 1910+9 m s _I . 
Furthermore, the sunspots embodying this regime were of a 
younger type. The second group, the scattered group, contained 
mostly old and perhaps recurrent sunspot groups. In the ve- 
locity, this group showed a large scatter with the mean speed 
of 1850 ms _1 . At that time we did not have a tool to follow 
a particular region on the Sun for any time of interest, there- 
fore we concluded that the existence of the above-mentioned 
regimes is compatible with the theory of dynamic disconnec- 
tion. Furthermore, the fact that all the young sunspots we ob- 
served depicted almost the same speed, led us to speculate that 
this behaviour is related to the assumption that the magnetic 
field emerges from a similar depth in the convection zone, in 
line with deep anchorin g. Following the appro x imatio n of the 
radial rotation profile in ICorbard & Thompsonl d2002l) . the ra- 
dius where the rotation corresponds to 1910+9 ms _1 is roughly 
0.946+0.008 R Q . 

In the past year, we have developed a tool for the selection 
of the active regions in the Michelson Doppler Imager (MDI; 
IScherrer et all [19951) magnetograms. From the previous studies 
we possess the datasets that cover all the MDI Dynamics cam- 
paigns in years 1996-2006, i.e., 502 days of high quality obser- 
vations. In each of these days, two full-disc 24-hours-averaged 
horizontal velocity maps, with an effective resolution of 60" 
sampled by 12 hours, were calculated. This allowed us to study 
the evolution of the active region dynamics recorded in our ve- 
locity datasets, in particular the phenomenon of the dynamic dis- 
connection from the magnetic roots. 

2. Data and Method 

In recent papers (e.g. ISvanda et al.L 120061) we introduced a 
method to measure the large-scale dynamics in the solar photo- 
sphere. This method is based on supergranular-structures track- 



ing in the full-disc, processed Dopplergrams, measured by MDI 
on-board the SOHO spacecraft. The application of this method 
allows one to compute the 24-hour averaged horizontal flow 
fields with resolution of 60" and noise level of 15 ms _1 . In 
the magnetised regions, this method makes it possible to mea- 
sure the apparent motion of supergranular-scale magnetic fea- 
tures (ISvanda et alll2009l) . 

The methodology consists of several steps. The main pur- 
pose of the pseudo-pipeline is to suppress disturbing effects, 
to remap the data onto a suitable coordinate system, to com- 
pute the horizontal vector-displacement field, and to convert the 
map of the displacements into the horizontal velocity vectors. 
The method and its subsequent validation are described in the 
great detail in lSvanda et alj d2006l) . Here we provide a very brief 
overview. 

The method processes 24-hour series of MDI full-disc 
Dopplergrams containing 1 440 frames. The one-day series first 
undergoes noise substraction and the removal of other distorting 
effects (Carrington rotation profile, p-modes). Then, the process- 
ing of averaged frames consists of two main steps. The first step 
involves the calculation of mean zonal velo cities using a very 
approximate local correlation tracking (LCT: iNovembeii 1 1 9861) 
algorithm. The obtained zonal velocities are fitted by a smooth 
function in a form of a> = cq + c\ sin 2 b + q sin 4 b, and the se- 
ries is tracked with the computed rotation profile. The tracking 
with the smooth velocity profile is done in order to minimise the 
average displacement of the supergranular structures caused by 
the solar differential rotation. Therefore the coefficients Co, c\, 
and C2 vary from series to series. In the second step, the LCT 
algorithm, with an enhanced sensitivity, is applied to obtain a 
low-noise-level displacement map. Finally, the differential rota- 
tion (obtained in the first step) is added to the vector velocity 
field obtained in the second step. Both steps can be divided into 
several sub-steps, which do not differ much in application. 

1 . The data series containing 96 averaged frames is tracked us- 
ing a selected rotation profile (Carrington rotation in the first 
step, differential rotation in the second step), and the frames 
are transformed into the Sanson-Flamsteed coordinate sys- 
tem to remove the geometrical distortion caused by the pro- 
jection onto the disc. The Sanson-Flamsteed (also known as 
"sinusoidal") pseudo-cylindrical projection conserves the ar- 
eas and is therefore suitable for the preparation of the data 
used by LCT. 

2. The tracked datacube undergoes k-co filtering with the cut- 
off velocity of 1500 ms _1 to suppress the noise coming 
mostly from the evolutionary changes of the supergranules. 

3. The LCT is applied: the lag between the correlated frames 
is 4 hours, the correlation window has a FWHM of 60", the 
measure of correlation is the sum of the absolute differences 
of the frames, and the nine-point method for calculation of 
the subpixel value of displacements is used. The calculated 
velocity field is averaged over the period of one day. The 
magnitudes of the vectors are corrected using formula 

^.corrected J ^ 1 3 ( ^ ?A,com P ute ^ \ _|15ms \ , ^ 

^;y,corrected / \ ^y,computed J \ / 

This formula results from the direct comparison between the 
synthetic data and the results of the displac ement measuring pr o- 
cedure applied to these synthetic data (see lSvanda et al.Ll2006l) . 

The processed dataset contains 1004 full-disc flow maps in 
502 days, when high-cadence MDI data were available and of a 
good quality. The velocity maps consist of many components on 
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Fig. 1. The Joy's law in our sample of active regions. Left - calculated for the sample of active regions, which emerged in CMD 
less than 60 degrees. Right - the same for active regions, that emerged in CMD less than 60 degrees and spent at least 4 days within 
CMD less than 60 degrees, so due to the averaging in time the measured tilts are less noisy. 



various spatial scales that cannot be reliably and unambiguously 
separated. On the largest scales, the differential rotation and the 
meridional circulation operate. We assume that these large-scale 
components vary slowly with the time, but their changes may 
influence the measurements of the proper motions of the active 
regions. Therefore, from each velocity map we subtracted the 
13-day running average to remove the systematic changes in the 
flows on the largest scales. We assume that other components 
within our resolution of 60" are the components of our inter- 
est. The removal of the long-term average will also suppress any 
possible systematic errors. 

Based on the daily bulletins coming from the Space 
Environment Center of the National Oceanic and Atmospheric 
Administration (SEC NOAA) we identified 564 labelled active 
regions in these maps. From which 522 were observed in more 
than one flow map. The theory predicts that the change in the 
dynamic regime should occur within a few (perhaps three) days 
after the magnetic regions has emerged in the photosphere. From 
the sample we therefore have to choose only those active regions 
that have emerged in the visible hemisphere within a reasonable 
distance from the central meridian, to avoid the influence of the 
edge effects. Only 194 active regions in the dataset were seen for 
the first time less than 60 degrees from the central meridian. This 
unfortunately means that the particular region could already be 
up to 1 day old at that time. Therefore, we included the mea- 
surements from the previous two days of the same region on the 
Sun into the analysis. From this set, we selected 72 bi-polar re- 
gions that survived at least 4 days within the central meridian 
distance (CMD) of less than 60 degrees. This drastic reduction 
of the sample was necessary in order to avoid the edge effects. 
Furthermore, to show that the disconnection-like behaviour is 
common among active regions showing significant changes in 
their dynamics over the lifespan, we chose another reduced sam- 
ple of active regions which possessed a standard deviation of the 
rotation speed of more than 10 m s _1 over the lifespan. We have 
18 such active regions in the sample. Note that the 10 ms _1 cri- 
terion is arbitrary. 

The 24-hour averaged velocity maps and corresponding 
magnetograms typically contain more than one active region 
at once. Thus, the selection of individual active regions is 
needed. This selection is based on the masks applied to the 
full-disc displacement maps and magnetograms. The masks for 
each active region were obtained in a semi-automatic way. 



Firstly, from the daily reports from the SEC NOAA, the loca- 
tion of the active region is determined and identified in MDI 
magnetograms. Then, for a wider area in the vicinity of the 
expected position of the active region, the structures of the 
magnetic field were roughly chosen by eye as the region-of- 
interest (closely corresponding with images from NASA's Solar 
Monitor, www.solarmonitor.org). The actual active region was 
identified in this region-of-interest by selecting pixels with a 
magnetic field strength above a threshold of 100 G. The 100-G 
region was then smoothed by 20", which created the final mask. 
This mask was then applied to the set of displacement maps and 
magnetograms corresponding to days when the particular active 
region was recorded in the SEC NOAA reports, and two days 
before. 

For each active region in the sample of 72, we computed the 
properties suitable for a detailed investigation. To study the evo- 
lution of flows in the active regions and the properties describ- 
ing these active regions, we calculated the large-scale flow field 
in the area covered by the magnetic field (with the threshold of 
100 G), and the total unsigned flux. As additional parameters, we 
computed the mean magnetic field intensities in both the leading 
and trailing polarities, the coordinates of the gravity centre of 
both polarities, and the mean proper motions of both polarities 
in the frame co-rotating with the mean rotation rate of the active 
region. All these properties are sampled by 12 hours and cover 
the whole interval, when the particular active region was located 
within a CMD of 60 degrees. 

We searched for any significant changes in the dynamics, 
which could be caused by the disconnection of the magnetic 
structures from their deeper roots. The expected change should 
reveal itself as a sudden deceleration of the proper motion of the 
active region, because the radial gradient of the mean plasma ro- 
tation is negative in near-surface subphotospheric layers. All the 
active regions in the sample underwent this investigation in or- 
der to obtain a homogeneous set of active regions properties to 
examine. 

3. Results 

3.1. Testing the sample 

The aim of this section is to show that the sample of active re- 
gions we have available displays the widely-accepted statistical 
properties of the magnetic field in the solar photosphere and is 
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therefore suitable for more detailed analysis of the proper dy- 
namics of active regions. We have selected two basic properties 
to test: the Joy's law a nd the polar i ty sep aration evolution. 

As first noted by lHale et all dl919l) . the bi-polar sunspot 
groups are tilted with respect to the zonal direction. The lati- 
tudes of the leading polarities are on average lower than the lat- 
itudes of the trailing ones. The so-called Joy's law shows that, 
statistically, the bi-polar regions are more tilted, if they are lo- 
cated at higher latitudes than the ones located at lower latitudes. 
The tilts are cons i dered essential by kinem atic dynamo mod- 
els (e.g. iLeightonl 1 19691 : I Wang et all 119911) and are often ex- 
plai ned as the action of the Cor iolis force on rising flux tube 
(e.g. lD'Silva & Choudhuriill993h . 

To test the Joy's law in our sample, we calculated the average 
magnetogram for each active region that emerged and remained 
in a CMD of less than 60 degrees. The tilt was measured from the 
centre-of-mass position of both polarities. The results are shown 
in Fig.Q] We clearly see that the tilts are mostly positive on the 
southern hemisphere and negative on the northern hemisphere, 
i.e., the leading parts of the bi-polar active regions are closer to 
the equator than the following ones. Joy's law itself is not nicely 
recove red. However, as pointed out in iKosovichev & Stenflol 
d2008l) . Joy's law holds only statistically and we cannot expect 
much better results than t hose presented in Fig, p] Our results 
are very similar to those in lKosovichev & Stenflol ( 20081) . Other 
effects, suc h as a depen d ence o n the phase of the solar cycle 
as noted by lUlrich et"ail (120021) may also take place. To sepa- 
rate such subtle effects the sample we have at our disposal is too 
sparse. 

Another property tested was the evolution of the separation 
distance between the leading and the following polarity in the 
spot g roup. Many works (e . g Ivan Driel-Gesztelvi & Petrovavl, 
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1 19901: iMoreno-Insertis et al.l 
pointed out that during the evolution, the separation between 
the leading and the following polarities grows. The active re- 
gion increases its size in the longitudinal direction. It is often 
interpreted as the result of the stretch caused by the differential 
rotation on the tilted region. The trailing polarity is at higher 
latitudes, therefore it senses a slower rotation than the leading 
part. We used our sample to measure this quantity. The separa- 
tion of the leading polarity from the following one is defined as 
the distance between the gravity centres of both polarities in the 
Carrington coordinate system. We constructed the histogram of 
the separation speeds based on the sample of 194 active regions. 
The results are shown in Fig. [2] We see that, on average, the sep- 
aration of polarities in the active regions is increasing with time, 
with the median speed of 0.3 heliographic degree per day, which 
corresponds to 45 ms~' on the solar equator. 

3.2. Acceleration or deceleration? 



As IRuzdjak et al.l d2004l) showed in a nice paper, based on the 
long-term observations of sunspots at Greenwich observatory 
covering more than a century, that rotation of sunspots usually 
slowed down during their evolution. I n fact, this observation 
was initially shown bv lTuominenl d 1 9621) . These papers primarily 
deal with the deceleration of sunspots in term of linear relation- 
ship with time. However, there is no principal reason to assume 
that in a highly dynamic layer, such a s is the convection zone , 
the deceleration should be monotonic. Siv araman et al. I (120031) . 
using Kodaikanal data, found an opposite behaviour, i.e. that the 
rotation rates of spot groups increase with their age. They sorted 
the sunspot groups according to their life spans from one to eight 
days, and found that the rotation of the groups on the first day 
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Fig. 2. The histogram of the polarity-separation speed for the 
sample of 194 active regions. It is evident that for most of the 
bi-polar active regions, the polarities are continuously separat- 
ing. The characteristic (median) speed is overplotted by vertical 
line and it is roughly 0.3 heliographic degree per day. 



of their life-cycle is the slowest. They then proceed to acceler- 
ate with age and reach the maximum velocity of rotation the day 
before dissolution. 

We studied the evolution of the proper velocities of the 
sunspots with respect to the averaged large-scale background, 
with a sampling of 12 hours. This sampling allowed for an 
improvement of the statistics compared to the simple average 
over the active region lifespan. Four different trends, which are 
present in the sample, are displayed in Fig. [3] The error-bars 
represent the statistical errors showing the l-cr scatter of the ve- 
locities over the whole active region. 

We observe that for two of these trends, i.e., continuous 
acceleration and variable behaviour, the change in the rotation 
speed is below 20 ms _1 . On the contrary, the other two cases 
display more significant changes. Table Q] contains the detailed 
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Fig. 3. Four different trends in the active region speed evolu- 
tion, which are present in our data sample. The regions are 
continuously decelerating (NOAA 8183), or continuously ac- 
celerating (NOAA 8038). Less clear trends are also noticed, 
such as the early acceleration followed by a rapid deceleration 
(NOAA 9408) and very variable behaviour (NOAA 8523). 
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Fig. 4. The division of the sample of the active regions into four dynamically different groups is based on the different shape of 
the trend of the active regions proper rotation speed with time. We divided 72 active regions in the group showing the continuous 
acceleration (upper left), continuous deceleration (upper right), acceleration followed by a sudden deceleration (bottom left), and 
the varying rotation speed (bottom right). With the thick line, the average trend among the displayed active regions is overplotted. 
Except for the third case, individual plots were not aligned. In the case with the possible dynamical disconnection imprint, we 
aligned the plots in time on the significant feature in the active regions lifetime, which is denoted by the one-to-two days lasting 
peak. 



results for all active regions in the sample. The sample is divided 
into four types based on the dynamic behaviour: The active re- 
gion is continuously accelerating (A), decelerating (D), shows 
the signs of the dynamic disconnection (DC), or behaves errati- 
cally (var). 

The tabulated results show that 33 % of the active regions 
in our sample first accelerate and then suddenly decelerate. In 
13 % of cases the active regions depict a continuous acceleration, 
while in 15 % of cases a continuous deceleration is observed. 
Finally, in the remaining 39 % a variable behaviour is present 
when the proper velocity does not change significantly or sys- 
tematically in the active region lifespan. The active regions ac- 
celerate mostly during the first few days of their lifespan. TableQ] 
also contains a number of other values describing the character- 
istic properties of the active regions. These include the measure 
of the average slope of the deceleration or acceleration in time 
for active regions falling in the corresponding category, the vari- 
ance of the rotation speed over the life-cycle of the active regions 
showing variable behaviour in their rotation speed. The measure- 
ments obtained for the active regions showing first acceleration 
and then deceleration can also be found in the same table, and 
are discussed in the following sections. 

Given the number of active regions in the sample, it is not 
feasible to demonstrate the exact trends in proper rotation evo- 
lution for each active region we studied. Instead, we associate 
these trends to the above-mentioned corresponding groups, with- 
out making explicitly clear which plot belongs to which active 
region. These plots can be seen in Fig. [4] and they demonstrate 



the overall trend in each selected group. Except for the third 
group, no alignment of the individual curves were done. The 
third group however contains a significant feature in the trend, 
a large peak lasting one to two days. The curves were therefore 
aligned in time using this feature. The thick solid line indicates 
the averages of all the trends contained in the panels. Using this 
average trend we are able to clearly distinguish different regimes 
for active regions showing continuous deceleration and for those 
representing the disconnecting group. Although it is difficult to 
clearly see the difference between the group with the assigned 
continuous acceleration and the variable behaviour, in the for- 
mer case one can see the overall increasing rotation speed with 
time when following individual trends. 

3.3. Possible evidence for disconnection 

From the sample, 18 active regions display a variance of speed 
evolution larger than 10 m s _1 , and from this 18, 78 % are of the 
type which display an initial acceleration followed by a sudden 
deceleration. We assume that these active regions are represen- 
tative of the sample which contain clear signatures of dynamic 
disconnection. 

We assume that disconnection is exhibited as a sudden de- 
crease in the proper speed, shortly before the maximum in area 
is reached. We further assume that the disconnection takes place 
at the instant of sudden decrease in speed. The upper part of the 
former flux-tube continues to rise and feeds more magnetic flux 
into the photosphere. Therefore the area of the magnetic field 
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Fig. 6. The example of the evolution of flows, active region area, the total flux, and the histogram of the magnetic field intensities. We 
see the overall acceleration in first two days (in this phase, the polarities separate rapidly), and sudden deceleration after this instant. 
The area of the active region continues to grow for another two days, the total flux in the active region for three days (however, 
usually the maximum in the area corresponds with the maximum in the total flux). In the bottom row of figures, in each column the 
histogram of the magnetic field strength in the region-of-interest in the particular part is displayed in by colours. The evolution of 
the dynamics seems very symmetrical, which is not true for the histograms of the magnetic field strength. The actual evolution of 
the flows in this particular active region is displayed in Fig. [7] 



in the photosphere keeps growing until the disconnected part is 
fully emerged, after which the flux starts to disperse and the ac- 
tive region area decreases. In Fig. [5]we display the speed at the 
assumed disconnection point (the disconnection speed, v^) and 



the speed in the following minimum (the relaxed speed, v le iaxX 
these numbers are also summarised in Table [TJ The thick lines 
emphasise the set of 14 active regions which display a signifi- 
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Fig. 7. The mosaic of the evolution of the flow field in the active region NOAA 8524. The background colours represent the large- 
scale magnetic field in the region, the grey arrows the the flows in the area. The differently coloured arrows in the magnetized areas 
represent the velocities that were used for the studies of the dynamics of this particular active regions. 



"O 

a> 

CD 
Q. 



O 
B 

o 
en 



1- 

I ■ 

I 



Disconnection speed 
Relaxed speed 



cant change in proper speed during their evolution, and which 
also exhibit the signs of disconnection. 

The typical drop of the speed during disconnection is 
~50 m s~'. This value corresponds well with the bi-modal distri- 
bution of the equatorial rotation rates in the presence of sunspots 
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Fig. 5. The jump in the speed, which can be caused by the dy- 
namic disconnection. For 23 active regions from the sample of 
72 we see that the speed in the pre-disconnection phase is by 
some 50 m s~' larger than in the post-disconnection phase. This 
value does not depend on the active region area. By the thick 
lines, active regions belonging to the subsample showing signif- 
icant change in the speed during the lifespan is emphasised. 



(with vel ocities cluster e d arou nd 1850ms -1 and 1910ms -1 ), as 
found bv lSvanda et al.1 ((2008). A quick check using the sunspot 
drawing archives showed that the "fast" group contained mostly 
the young and growing regions (which corresponds to the pre- 
disconnection phase in the current study), while the second 
"scattered" group contained mostly old and dispersing regions. 
In the current study this corresponds to the post-disconnection 
phase. 



Although our analysis is based on the mean rotation rate con- 
sidered for the whole active region, the signatures of the dis- 
connection should be detectable in both polarities. However, it 
is difficult to measure this phenomenon, because the motions 
of polarities in the co-rotating frame are an order of magnitude 
faster. The dynamics in the co-rotating frame (as displayed e.g. 
in Fig.|6]l is usually very symmetrical in both polarities. Taking 
a careful look at this issue we found that in some 70 % of ac- 
tive regions in the sample showing clear signatures of the dis- 
connection, the disconnection features can be detected in both 
polarities, while in 30 % it is clear only in one polarity. For a 
few cases, the time when disconnection occurs in opposite po- 
larities differs by approximately 1 day. This issue is probably 
related to the configuration and asymmetry of the magnetic field 
in the particular active region. The disconnection is proposed for 
flux-tubes that form sunspots or pores, this does not apply to the 
magnetic field forming the plage. Unfortunately, as our method 
allows the measurement of large-scale features, we cannot ad- 
dress this issue properly, and it can possibly affect our results. 
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3.4. Disconnection depth 

We observe a systematic shift f] ag between the time of the dy- 
namic regime change and the time of the maximum in the area 
(see e.g. Fig.|6]l. We interpret this lag as the time which the flux 
tube forming the magnetic island rises from its parental mag- 
netic structure to the surface after the disconnection. After the 
disconnection takes place, no more magnetic field is fed into the 
magnetic island and after its entire emergence the photospheric 
magnetic field starts to diminish. In all but two cases this lag 
is positive, therefore fulfilling the above mentioned assumption 
(see TableQ}. 

The average lag is 1.6 ± 1.1 days. The flux tube rises in 
the convection zone with a speed, Vri se , which is unknown and 
must be determ ined for each cas e, e.g., by numerical simulations 
(such as in e.g. iFan etallfl993t) . We may roughly estimate this 
rising speed by various characteristic speeds in the convection 
zone, to put limits on the disconnection depth. Various speeds in 
the upper 70 Mm of the convection zone are displayed in Fig. [9] 

The most reasonable estimate of the rising speed is the 
Alfven speed, ca, 



4000 



B 



(2) 



because the perturb ations in the m agnetic field usually propagate 
with this speed (e.g. lParkerll 19751) . In the above equation B is the 
intensity of the magnetic field, is the permeability, and p is the 
plasma density. Based o n the solar plasma parameters contai ned 
in the standard Model S dChristensen-Dalsgaard et aU 1 19961) at- 
mosphere, we may estimate the rising time fn Se of the flux tube 
from the given depth in the convection zone /? start onto the sur- 
face: 



V^p(^) 

B 



dr. 



(3) 



Assuming a constant 1-kG flux-tube emerging with the 
Alfven speed, the theoretical fi ag is 180 days if the flux-tube 
emerges from the bottom of the convection zone, and 7 days 
when emerging from the subsurface shear at 0.95 R Q . Even the 
10 5 -G field would rise from the base of the convection zone for 
approximately one month. Although these numbers represent a 
very rough estimate, they show that a measured time-lag slightly 
more than 1 day corresponds to the disconnection location in the 
shallow subsurface layers. 

In Fig. [8] we show the theoretical rising time depending on 
the initial depth and the strength of the flux tube that is con- 
sidered constant during the rising process. The measured lag is 
displayed in contours. We have to consider our theoretical results 
as an estimate, because the calculation does not include the flux- 
tube expansion, which naturally influences the rising speed. The 
precise modelling and the reproduction of the situation would 
be very difficult, because the initial properties of the flux tube 
deep in the convection zone are unknown. We believe that our 
estimate shows the essence of the problem. Inverting ([3} with 
frfse ~ flag we can estimate the disconnection depth for the par- 
ticular active region. 

The Alfven speed gives the smallest disconnection depth 
estimate. Buoyant ascent of the flux tube requires the super- 
equipartition field, e.g., ca > v c , where v c is a convection speed. 
Using the convection speed as the rising speed after the assumed 
disconnection in our study therefore provides the largest discon- 
nection depth estimate. It is clear that the presence of the strong 
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Fig. 8. The theoretical time-shift between the disconnection time 
and the time when the maximum in area is reached, based on 
the assumption that flux-tube rises with the Alfven speed as the 
function of depth and flux-tube strength. The measured time- 
shift including the error is displayed in contours. This very rough 
estimate shows that the disconnection cannot occur deeper than 
some 40 Mm for 4-kG flux-tube to reproduce the observed be- 
haviour. The Alfven speed computa tion is based on Model S of 
IChristensen-Dalsgaard et af] dl996l) . The real rising speed of an 
expanding flux-tube should not differ from the Alfven speed in 
the order of magnitude. 



magnetic field alters convection and therefore probably changes 
the convection speed. However, the main purpose of this calcu- 
lation is to bound the disconnection depth, therefore the use of 
the model convection speed without the influence of the mag- 
netic field is justified. Again, without proper modelling it is not 
possible to go beyond prese nted estima tes. 

The mixing-length (e.g. lStix[ [l989) convection speed in the 
Sun can be computed from a known solar model (e.g., Model S 
in our case) using the following formula: 



v,- = a/- 



(4) 



where m(r) is the total mass in the sphere with radius r, a is the 
mixing-length parameter, H p is a pressure scale height, and A* 
is a convection parameter, 



A\r) = 



1 dlnp(r) dlnp(r) 



dlnr 



y{r) d In r 



(5) 



where y, p, and p represent the adiabatic exponent, pressure, and 
density. The rise of the flux-tube with a convection speed implies 
that the strength of the field is the equipartition one. 

In Table Q] we provide the range of the disconnection depths 
for each active region studied. We observed that the disconnec- 
tion occurs within a few tens of Mm in depth. The maximum 
intensity of the magnetic field used for the calculation of the 
Alfven speed is determined from the MDI magnetograms. We 
need to keep in mind that the actual strength of the original flux 
tube (also the rising time) may be different from the estimated 
one. 



4. Conclusions 

We conclude that the observed dynamical behaviour of a signif- 
icant fraction of active regions in the homogeneous sample in 
the past solar cycle is compatible with a dynamic disconnection 
of the bi-polar active regions from their parental magnetic roots. 
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In our sample, 33 % of active regions display the signs of the 
disconnection. In the reduced sample with significant changes in 
the dynamics, 78 % show the desired behaviour. The disconnect- 
ing active regions display a lag between the disconnection time 
and the instant when the maximum in area is established. If we 
assume that this lag is caused by the continuous rise of the dis- 
connected magnetic field up to the photosphere, we may estimate 
the depth where disconnection takes place. We estimate that dis- 
connection takes place within few tens of Mm below the solar 
surface, not in disagreemen t with 1-D simulations conducted by 
ISchiissler & Rempell d2005l) . 

Th e numerical simulations done by ISchiissler & Rempel 
d2005l) predict that the disconnection should happen for the sig- 
nificant fraction of bi-polar active regions, if not for all of them. 
However, our results show the clear attributes of the disconnec- 
tion in only one third of the studied cases. We speculate that in 
other cases, disconnection does not occur during the observed 
span (e.g. before emergence into the photosphere in the case of 
continuously decelerating active regions) or that the influence 
of the proper motions of the magnetic features by the turbulent 
background is too high. Therefore we cannot reliably separate 
the attributes of the disconnection from other motions. 

The mechanism of the dynamic disconnection is proposed 
for strong flux-tubes that form sunspots or pores. With our in- 
vestigation method we unfortunately do not have sufficient reso- 
lution to distinguish between the motion of magnetic structures 
in, e.g. sunspots, from the ones in the plage region, which prob- 
ably embody different dynamic properties. 

Our present study does not provide a definitive answer on 
the question of whether the origin of magnetic activity is in the 
deep or shallow dynamo action. Both approaches require some 
form of disconnection of the surface magnetic field to allow the 
dispersed flux to migrate towa rds the solar poles. The dee p dy- 
namo based model assumed bv lSchussler & Rempell ( 120051) pro- 
vides predictions for particular active regions dynamics, which 
we have tested. Although the results of this test favour the model 
assumed, they are not robust enough to rule out other interpreta- 
tions. The phenomenon of dynamic disconnection requires fur- 
ther investigation using high resolution methods and by local 
helioseismology, which proves its power in investigating deep 
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Fig. 9. The trends in various speeds in the upper part of the solar 
convection zone. The Alfven speed and the convection speed are 
used to estimate the limits in depths, where the disconnection 
takes place. The computation of the speed trends is based on 
model S. 



structure of sunspots (e.g. lMoradi et al.l , r2009t) . Clear, deep mag- 
netic field signatures would rule out a group of models assuming 
a shallow dynamo and provide other tests for rising flux-tubes 
models. 

Acknowledgements. The authors were supported by the Grant Agency of 
Academy of Sciences of the Czech Republic under grant IAA30030808. 
The Astronomical Institute of ASCR is working on the Research project 
AV0Z 10030501 (Academy of Sciences of CR), the Astronomical Institute of 
Charles University on the Research program MSM002 1620860 (Ministry of 
Education of CR). SOHO is a project of international cooperation between ESA 
and NASA. The authors thank the referee Roger K. Ulrich, whose comments and 
suggestions significantly improved the quality of the study and its presentation, 
and also Hamed Moradi for correcting the English. 



References 

Balthasar, H., Vazquez, M., & Woehl, H. 1986, A&A, 155, 87 
Bumba, V. & Howard, R. 1965, ApJ, 141, 1502 

Christensen-Dalsgaard, J., Dappen, W., Ajukov, S. V., et al. 1996, Science, 272, 
1286 

Corbard, T. & Thompson, M. J. 2002, Sol. Phys., 205, 211 

D'Silva, S. & Choudhuri, A. R. 1993, A&A, 272, 621 

D'Silva, S. & Howard, R. F. 1994, Sol. Phys., 151, 213 

Fan, Y, Fisher, G. H., & Deluca, E. E. 1993, ApJ, 405, 390 

Fan, Y., Fisher, G. H., & McClymont, A. N. 1994, ApJ, 436, 907 

Hale, G. E., Ellerman, F, Nicholson, S. B., & Joy, A. H. 1919, ApJ, 49, 153 

Herdiwijaya, D. 2002, in 8th Asian-Pacific Regional Meeting, Volume II, ed. 

S. Ikeuchi, J. Hearnshaw, & T. Hanawa, 427^127 
Kosovichev, A. G. & Stenflo, J. O. 2008, ApJ, 688, LI 15 
Leighton, R. B. 1969, ApJ, 156, 1 

Moradi, H., Hanasoge, S. M., & Cally, P. S. 2009, ApJ, 690, L72 
Moreno-Insertis, F, Caligari, P., & Schuessler, M. 1994, Sol. Phys., 153, 449 
November, L. J. 1986, Appl. Opt., 25, 392 
Parker, E. N. 1975, ApJ, 198, 205 

Ruzdjak, D., Ruzdjak, V, Brajsa, R., & Wohl, H. 2004, Sol. Phys., 221, 225 
Schatten, K. H. 2009, Sol. Phys., 255, 3 

Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 
Schiissler, M. & Rempel, M. 2005, A&A, 441, 337 

Sivaraman, K. R., Sivaraman, H., Gupta, S. S., & Howard, R. F. 2003, Sol. Phys., 
214, 65 

Sobotka, M. & Roudier, T. 2007, A&A, 472, 277 

Stix, M. 1989, The Sun. An Introduction (Springer- Verlag Berlin Heidelberg 
New York) 

Svanda, M., Klvana, M., & Sobotka, M. 2006, A&A, 458, 301 

Svanda, M., Klvana, M., Sobotka, M., & Bumba, V. 2008, A&A, 477, 285 

Svanda, M., Klvana, M., Sobotka, M., Kosovichev, A. G., & Duvall, Jr., T. L. 

2009, New Astronomy, 14, 429 
Tuominen, J. 1962, Zeitschrift fur Astrophysik, 55, 1 10 
Ulrich, R. K, Evans, S., Boyden, J. E., & Webster, L. 2002, ApJS, 139, 259 
van Driel-Gesztelyi, L. & Petrovay, K. 1990, Sol. Phys., 126, 285 
Wang, Y.-M., Nash, A. G., & Sheeley, Jr., N. R. 1989, Science, 245, 712 
Wang, Y.-M., Sheeley, Jr., N. R., & Nash, A. G. 1991, ApJ, 383, 431 



10 M. Svanda et al.: Large-scale horizontal flows in the solar photosphere V: A disconnection of bi-polar sunspot groups 



Table 1. The studied properties of 72 active regions in our sample. Active regions of type DC show signatures of dynamic disconnection. 



NOAA 


From 


To 


Type 


Max. intensity 


Vdis 


Vrelax 


Slope 


Variance 


flag 


Disc, depth 










Gauss 


ms _1 


ms" 1 


m s _1 day -1 




days 


Mm 


0651 


2004/07/13 


2004/07/22 


var 


700 








5.7 






0654 


2004/07/24 


2004/08/03 


DC 


100 


16(3) 


-7(2) 






-1.0 


(5.0-26.7) 


0782 


2005/06/27 


2005/07/05 


DC 


1300 


14(6) 


-14(4) 






-1.5 


(19.0-34.5) 


0783 


2005/06/29 


2005/07/09 


DC 


1300 


31(9) 


-6(2) 






-2.0 


(21.9-41.5) 


0785 


2005/06/28 


2005/07/06 


var 


800 








6.5 






0791 


2005/07/22 


2005/07/31 


A 


1400 






0.7 








0793 


2005/07/27 


2005/08/04 


DC 


700 


29(9) 


0(0) 






-1.0 


(11.7-26.7) 


7968 


1996/06/01 


1996/06/11 


var 


200 








3.7 






7978 


1996/07/02 


1996/07/12 


var 


700 








4.3 






8032 


1997/04/14 


1997/04/19 


var 


500 








4.7 






8038 


1997/05/06 


1997/05/16 


A 


700 






1.9 








8040 


1997/05/15 


1997/05/25 


DC 


800 


18(7) 


-8(2) 






-2.0 


(17.2-41.5) 


8048 


1997/05/29 


1997/06/09 


var 


1000 








15.5 






8050 


1997/06/08 


1997/06/18 


var 


300 








4.3 






8052 


1997/06/10 


1997/06/19 


A 


400 






0.3 








8059 


1997/07/02 


1997/07/10 


D 


200 






-1.3 








8131 


1998/01/10 


1998/01/18 


D 


1300 






-1.7 








8141 


1998/01/18 


1998/01/26 


var 


300 








3.6 






8142 


1998/01/20 


1998/01/28 


var 


800 








4.2 






8145 


1998/01/24 


1998/02/02 


A 


300 






0.1 








8153 


1998/02/06 


1998/02/13 


D 


400 






-0.6 








8156 


1998/02/11 


1998/02/21 


D 


1200 






-2.2 








8160 


1998/02/13 


1998/02/20 


var 


200 








4.8 






8161 


1998/02/16 


1998/02/21 


A 


300 






0.7 








8164 


1998/02/21 


1998/03/01 


A 


1000 






-0.0 








8167 


1998/02/22 


1998/03/02 


DC 


800 


19(6) 


-2(2) 






-0.5 


(9.0-17.4) 


8174 


1998/03/05 


1998/03/15 


var 


1000 








4.0 






8175 


1998/03/07 


1998/03/16 


var 


200 








6.7 






8180 


1998/03/11 


1998/03/21 


D 


800 






-1.1 








8183 


1998/03/18 


1998/03/28 


D 


100 






-2.4 








8184 


1998/03/19 


1998/03/26 


var 


800 








3.6 






8188 


1998/03/25 


1998/04/02 


DC 


200 


42(11) 


-24(4) 






-4.0 


(12.4-64.9) 


8191 


1998/03/30 


1998/04/07 


var 


500 








6.6 






8193 


1998/04/02 


1998/04/09 


var 


1800 








4.0 






8507 


1999/04/03 


1999/04/12 


var 


200 








5.4 






8509 


1999/04/05 


1999/04/14 


var 


500 








8.0 






8511 


1999/04/06 


1999/04/16 


var 


1300 








4.2 






8512 


1999/04/07 


1999/04/17 


var 


800 








6.4 






8513 


1999/04/08 


1999/04/17 


var 


800 








4.1 






8515 


1999/04/11 


1999/04/19 


D 


800 






-1.2 








8523 


1999/04/23 


1999/05/02 


DC 


300 


3(2) 


-8(3) 






-3.0 


(13.1-53.9) 


8524 


1999/04/23 


1999/05/02 


DC 


1200 


42(19) 


2(5) 






-2.0 


(21.0-41.5) 


8526 


1999/04/27 


1999/05/07 


var 


200 








3.2 






8530 


1999/05/01 


1999/05/12 


D 


500 






-1.3 








8536 


1999/05/05 


1999/05/14 


var 


1100 








5.0 






9024 


2000/05/28 


2000/06/07 


DC 


900 


28(10) 


-4(2) 






0.0 


(0.0-0.0) 


9028 


2000/05/31 


2000/06/08 


var 


400 








3.7 






9030 


2000/06/03 


2000/06/13 


A 


500 






-0.3 








9032 


2000/06/03 


2000/06/13 


DC 


600 


12(2) 


-5(3) 






-4.5 


(22.3-70.0) 


9035 


2000/06/07 


2000/06/15 


DC 


500 


10(5) 


-5(3) 






-0.5 


(7.4-17.4) 


9037 


2000/06/09 


2000/06/19 


DC 


600 


10(3) 


-8(2) 






-1.5 


(13.1,34.5) 


9050 


2000/06/15 


2000/06/23 


DC 


500 


16(10) 


-9(3) 






-1.0 


(10.0-26.7) 


9059 


2000/06/21 


2000/06/29 


var 


500 








5.9 






9063 


2000/06/25 


2000/07/05 


var 


600 








3.0 






9365 


2001/03/01 


2001/03/06 


DC 


500 


46(10) 


0(0) 






0.0 


(0.0-0.0) 


9400 


2001/03/23 


2001/03/30 


A 


400 






5.3 








9401 


2001/03/25 


2001/04/03 


D 


1300 






-8.1 








9408 


2001/03/27 


2001/04/03 


DC 


1600 


38(7) 


-4(4) 






-1.5 


(21.0-34.5) 


9424 


2001/04/08 


2001/04/16 


DC 


700 


33(14) 


-2(3) 






-1.0 


(11.7-26.7) 


9426 


2001/04/10 


2001/04/18 


var 


700 








7.1 
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Table 1. continued. 



INUAA 


From 




lype 


Max. intensity 


Vdis 


^relax 


L' 1 ..... . 

Mope 


Variance 


flag 


Disc, depth 










Gauss 


_i 

ms 


_i 

ms 


m s day 


_i 

ms 


days 


Mm 


9431 


2001/04/16 


2001/04/26 


DC 


700 


23(5) 


-18(5) 






-1.0 


(H.7-26.7) 


9432 


2001/04/16 


2001/04/24 


DC 


900 


19(8) 


-10(2) 






-1.5 


(15.9-34.5) 


9441 


2001/04/26 


2001/05/05 


DC 


1200 


21(11) 


-3(1) 






-2.5 


(23.5-47.9) 


9442 


2001/04/26 


2001/05/06 


D 


400 






4.7 








9910 


2002/04/17 


2002/04/24 


DC 


1200 


25(15) 


-3(2) 






-2.5 


(23.5-47.9) 


QQ1 S 

771J 


9009/n/L/99 
ZUUZ/LHyZZ 


9009/04/90 


VST 


ljUU 








j.j 






9922 


2002/04/22 


2002/05/01 


A 


200 






-2.0 








9932 


2002/05/01 


2002/05/11 


var 


500 








4.0 






9935 


2002/05/01 


2002/05/10 


D 


500 






0.4 








9936 


2002/05/01 


2002/05/10 


DC 


400 


13(4) 


-6(3) 






-1.0 


(9.0-26.7) 


9950 


2002/05/10 


2002/05/19 


var 


700 








5.8 






9970 


2002/05/26 


2002/06/02 


DC 


1500 


0(3) 


-10(3) 






-0.5 


(12.0-17.4) 



